Exact Analysis of ESR Shift in the Spin-1/2 Heisenberg Antiferromagnetic Chain 
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A systematic perturbation theory is developed for the ESR shift and is applied to the spin- 
1/2 Heisenberg antiferromagnetic chain with a general anisotropic exchange interaction. Using the 
Bethe ansatz technique, the resonance shift is obtained exactly for the whole range of temperature 
and magnetic field, in the first order of the anisotropy. The obtained g-shift strongly depends on 
magnetic fields at low temperature, showing a significant deviation from the previous classical result. 
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Many of the interesting phenomena and experimental 
measurements in strongly correlated quantum systems 
are related to dynamics of the system. Despite the con- 
siderable progress in many-body theory, available exact 
results on quantum dynamics in many-body systems are 
still rather scarce. Indeed, even a systematic framework 
for approximate calculations is often not well established. 
Electron spin resonance (ESR) is a typical example of 
such dynamical phenomena . 

In ESR, a static magnetic field H and an electromag- 
netic wave with the polarization perpendicular to H are 
applied to a magnetic system, and then the intensity of 
the radiation absorption /(w) is measured as a function of 
frequency to of the electromagnetic wave. Within the lin- 
ear response theory, I{uj) is written as /(w) cx ujx'!rcr{'^)^ 
where Xaa (^) is the imaginary part of the dynamical sus- 
ceptibility Xcrcr(w) at zero momentum: 



[5"(0,5"(0)])e*"* dt. 



(1) 



Here N is the total number of spins, 5'°' is the total 
spin operator {S" — ^2^=1 ^j) ^ ^hc direction 
of the polarization. We take the 2;-axis as the direction 
of the magnetic field H, and for simplicity we assume 
a G {x, y}. Let us consider a system with Hamiltonian 



SU{2) 



Hz + Ti.' , 



(2) 



where Ti.su{2} is an >S'[/(2)-symmetric term and Tiz = 
—HS^ is the Zeeman term, and Ji' is an anisotropic term 
assumed to be a small perturbation. In the absence of 7i', 
the energy absorption occurs only at the paramagnetic 
resonance frequency uj — H . It is only the anisotropic 
perturbation 7i', although it may be small, that causes 
the change in the spectrum such as a shift and a broad- 
ening of the resonance peak. 

Several frameworks of perturbation theory in 7i', typ- 
ified by the Kubo-Tomita theory 0], have been devel- 
oped. However, the questionable assumptions and tech- 
nical difficulties limited their success in systems with 
strong quantum fluctuations. More recently, field theo- 
retical 0] and numerical Q approaches avoiding ambigu- 
ous assumptions have been provided and brought about a 



renewed interest in ESR theory. Unfortunately, however, 
the former is applicable to ID systems only in low-energy 
regimes, and the latter is applied only to small systems. 
Hence rather much remains open in ESR theory. 

This is evident in theory on the shift in the Heisenberg 
antiferromagnetic chain, which is one of the most fun- 
damental problems in ESR. The ESR shift 5u) is usually 
defined as the deviation of the resonance frequency from 
the paramagnetic resonance frequency oj = H . Nagata- 
Tazuke theory ^ published more than 30 years ago has 
been the only theory for arbitrary temperature, and still 
remains the standard. However, it is based on the clas- 
sical spin approximation and does not incorporate the 
quantum fluctuation. In this letter, we utilize the in- 
tegrability of the the spin-1/2 Heisenberg antiferromag- 
netic chain to go beyond the limitations of the previous 
approaches. As a result, we obtain the ESR shift for 
general anisotropic interactions at arbitrary temperature 
and magnetic field, fully including the quantum effects. 

Nagata-Tazuke theory 5] is based on the fundamental 
formula for the ESR shift, first proposed in ref. Its 
simplest form reads: 



5lj = — 



{[AS-]) 

2{S'-) ' 



(3) 



where (• • • ) is the thermal expectation value with respect 
to H, and A = [H', 5+]. In ref. [3, this formula is further 
evaluated for the ID classical Heisenberg antiferromagnet 
in the lowest order of H, utilizing Fisher's solution 7] of 
the classical Heisenberg chain. Owing to the reliance on 
the classical limit, we cannot expect their result to apply 
to the strongly quantum system. On the other hand, as 
some nontrivial assumptions are involved in the original 
derivation of the formula Q, in literatures there seems 
to be a confusion over the validity of the starting point, 
eq. Q itself. 

Here, we clarify the validity and limitation of the for- 
mula To discuss the resonance shift unambiguously, 
in this letter, we define the shift 6ujaa for the polarization 
a in terms of the dynamical susceptibility as 



H. 



(4) 
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Note that the integration interval is taken to be [0, oo] 
instead of [— cx),oo], because Xa^i^) is an odd function 
of a;. By using a sum rule the numerator of eq. Q 
is explicitly written as —tt{[['H,S'^],S''])/2. In contrast, 
the integration in the denominator cannot be performed 
analytically. Hence we shall evaluate (0J perturbatively. 

The computation in the first order can be carried out 
without resorting to any other approximation, by uti- 
lizing the exact Heisenberg equation of motions such as 
S+ = -iHS+ + iA^. After straightforward but some- 
what tedious calculations, we obtain the first order per- 
turbative contribution to eq. Q, exactly as in eq. 
Namely, eq. (jSJ is generally exact in the first order of the 
anisotropic perturbation Ti.' . Note that, however, ei 
breaks down in the second and higher orders in Ti' 
Restricting eq. to the first order in Ti.' , we obtain the 
simplified expression 



{[A,S-])o 



(5) 



where (•••)o is the thermal expectation value with re- 
spect to the unperturbed term Ho = Hsu (2) + T^z- 

As there is no contribution to ^ from satellite 
peaks , it exactly gives the shift of the main paramag- 
netic resonance in the first order. As we will see, for the 
exchange anisotropy (e.g. dipolar interaction), the first 
order contribution is generally non-vanishing. Thus the 
precise evaluation of eq. © is of great importance. 

Now, let us consider the Heisenberg chain 



N 



N 



T^o = 'Hsu (2) +'H.z ^ j'^Sj ■ Sj+i - H^^Sj, (6) 

with a small anisotropic interaction as a perturbation 
Ti.' . As Ti', we consider the anisotropic exchange in- 
teraction. Since the antisymmetric exchange interaction 
(Dzyaloshinskii-Moriya interaction) does not contribute 
in the first order, we only consider the symmetric (with 
respect to an exchange of neighboring spins) exchange 
interaction. The dipolar interaction may also be approx- 
imated by a suitable exchange anisotropy in the nearest 
neighbor interaction. Assuming the chain is uniform, the 
anisotropic exchange interaction can always be diagonal- 
ized as 



j p,q£{a,b,c} 



T/ CP qq 



(7) 



where a,b,c are the principal axes of the anisotropy. 
Let the direction of the static magnetic field H be 
(q!,/3, 7) = (sin cos 0, sin sin 0, cos 0) in the (a, 6, c)- 
coordinate, where 6 is zenith angle and (/> is azimuth an- 
gle. As introduced previously, the (x, y, z)-coordinate is 
defined so that H is applied in the z-axis. The resonance 
shift (El is now evaluated as 



Here the direction-dependent factor f{9, (/>) and the field- 
and temperature-dependent factor Y{T,H) are respec- 
tively given by 

fie, ^) - j:(I - 3a') + 4(1 - 3/3^) + J'^il - 37'), (9) 



Y{T,H) 



QX QX \ 



(5|)o 



(10) 



In general, eqs. are valid for arbitrary spin, when 

we consider Ti' as a perturbation. It would be inter- 
esting to evaluate them numerically, for example by a 
Quantum Monte Carlo simulation. However, this could 
be a challenging problem especially in the weak field 
regime; there will be a cancellation of significant digits 
in taking the difference between the very close quantities 
and (5/5^1). 
Hereafter we shall concentrate on the 5=1/2 case. 
The total Hamiltonian (0) is no longer integrable even 
with S = 1/2, due to the existence of the anisotropy 
and magnetic fields. The first-order shift 1)10(1 . how- 
ever, can be exactly calculated by the Bethe ansatz tech- 
nique [iq| utilizing the integrability of the zeroth order 
Hamiltonian 10. After long calculations, we obtain 

Y{T,H) = ]^~^j>^\n{l + 7^{x + i))dx, (11) 

where the contour T encloses the real axis counterclock- 
wise: for instance T can be the rectangular contour whose 
corners are —00 — z, +00 — i, +00 + i, and —00 -I- i. 

The unknown function ?/(x) is determined by the fol- 
lowing non-linear integral equation: 

\nr]{x) = ^^ai{x) - ^ 

a2{x-y-i) \n{l + ri{y + i)) dy, (12) 

where a„(a;) = n/(7r(x^ + n^)). 

Although the result itTT)) is hmited only to the S* = 1/2 
Heisenberg antiferromagnctic chain, it gives the reso- 
nance shift H10|) without the difficulty due to the can- 
cellation of significant digits. As a consequence, the res- 
onance shift can be evaluated numerically with a very 
high accuracy for arbitrary fields and temperatures. 

In experimental studies, it is customary to discuss the 
shift in terms of the effective g-factor 



5cff = 5oo 



9o 



fie,4>)Y{T,H) 
H 



(13) 



6lo = f{0,^)xY{T,H). 



(8) 



and the g-shift Ag = gcs — goo ■ Here Su> is the frequency 
shift we have obtained with setting ^bQoo — 1, and 1700 
is the g-factor appearing in the Zeeman term. When 
Jp <C J, H, we can identify goo with the observed g-factor 
in the high temperature limit. 

In Fig.^ the temperature dependence of JY{T, H) / H, 
which is proportional to the (7-shift Ag, is depicted for 
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FIG. 1: The temperature dependence of JY{T, H) / H for var- 
ious H. For comparison, the result from the classical approxi- 
mation is also depicted. Inset: details of the low-temperature. 



FIG. 2: Comparison of the exact result with the asymptotic 
behavior H16^ (we set C = 3/4). A crossover from a regime 
dominated by logarithmic corrections to a regime dominated 
by algebraic corrections occurs aX T ^ H. 



various magnetic fields. We find that Ag clearly de- 
pends on the magnetic fields H especially in the low- 
temperature regime, and converges to a certain value at 
T = for any finite H . One also finds that Aji decreases 
with increasing H. In the limit H oo, Ag vanishes 
at any temperature. In contrast, for a weak magnetic 
field <C 1, Ag logarithmically increases with decreas- 
ing temperature and approaches infinity when — > 0. 

In Fig. m the present results are compared to the pre- 
vious classical result 0. It was obtained by expand- 
ing eq. l(Tn|l in the lowest order of H and then applying 
Fisher's exact solutions to the classical Heisenberg chain 
at ff = 0, as 



Yci{T, H) 



HK r 2 



5J 



1 



u{K)lK 2K 



u{K) = cothX - 



JS{S 



K = - 



1) 



T 



(14) 



where we shall set S — 1/2 in the present case. Obvi- 
ously the g-factor for the classical approach diverges with 
decreasing T, which contrasts to our results. In addition, 
the classical result does not describe the i7-dependence 
as it is obtained in the lowest order expansion in H . 

Comparing the integral equations Hll(l and 1)121) for 
Y{T,H) and that for the magnetization (S'^)o .10, ll| 
in the S* = 1/2 Heisenberg chain ((HJ, we find a remark- 
able identity relating them: 



Y{T,H) = ^ I {SI 



dJ. 



(15) 



This is useful in obtaining the low-temperature asymp- 
totics of the shift. Because the low-temperature asymp- 
totics of {S^)a is derived by an effective field theory [l2| . 



we obtain, for < T < J, 



y(r, H) 



H 
1^ 



+ c}, c 



InC 



1 

AC 



2j27r 



In 



V \/2T 



(16) 



Here is the Euler constant, and the integration con- 
stant C ~ 3/4 is determined by the fitting to the exact 
result eq. Hll|) , as in Fig. [3 

Note that the leading log in H16|) agrees exactly to what 
was obtained in ref. |3| via a field theoretical approach. 
However, the "leading log" is only dominant in an ex- 
tremely low temperature and magnetic field. For a real- 
istic parameter regime, the subleading (in J/T or J / H) 
corrections determined in the present study is actually 
essential in a quantitative discussion of the shift. 

Furthermore, our analysis determines the running cou- 
pling constant with respect to a marginal operator in 
ref. [sj precisely. Following the logic of ref. .3,j , this allows 
a prediction of the linewidth n including the subleading 
corrections for 7^ 0, = = 0: 



7^3 



1 



ln£ 



C 



} , (17) 



where e = 4 (e = 2) when the applied field H is parallel 
(perpendicular) to the c-axis. The subleading corrections 
should be also essential in a quantitative prediction for 
realistic values of temperature and magnetic field. 

Finally, we discuss experimental data in the light of 
our results. In Fig. |31 we analyze the g-shift data on 
LiCuV04, which is considered to be an 5* = 1/2 Heisen- 
berg antiferromagnetic chain. The data are taken from 
ref. ^3 (powder sample) and ref. (single crystal). 
We take the estimate J = 44 [K] from ref. 0, 01 for 
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the intrachain isotropic exchange interaction Then 
H/J = 0.027 and H/J = 0.010, respectively, for ref. [H 
and ref. Q. For these parameters, in Fig. |3|we plot the 
theoretical result Y{T,H) from eq. The difference 

in the magnetic field H is neghgible in the temperature 
regime in which we discuss the experimental data. 

The principal axes of the g-tensor, i.e. {g°',g^,g'^) and 
those of the anisotropic exchange tensor (O are thought 
to coincide with the crystallographic axes (a, 5, c) [ij, 
Q|. Therefore, we use the bare factor g^^''^ and the 
anisotropic exchange J'^^i, ^ for each crystallographic axis 
as fitting parameters. Here we note that we can set one 
of ^ to zero, by renormalizing the isotropic exchange 
J. The normalized (7-shift corresponding to Y{T, H) was 
then obtained for each set of experimentally observed 
effective g- factor, using eqs. ^ and ((Tni) . 

In ref. the anisotropy in a5-plane was neglected 
in analyzing the powder sample. For this set of data, we 
obtain g^ = 2.22, = g'^ = 2.09 as a result of fitting in 
Fig. 13 The best fit value of the anisotropic exchange was 

0.15 J fori? II cand O.llJ for iJ _L c, where 

we set J'a ~ Jb ^ 0. For the data from ref. (IJI, our fitting 
yields = 2.31,5^ = 2.07, = 2.09. The set of the 
anisotropic parameters ~ — 0.085J, ^ 9 x 10^"^ J 
and = well reproduces experimental data for the 
all directions H \\ a,b,c. These estimates are somewhat 
different from those in ref. ^3|- In any case, the observed 
small anisotropy in the a6-plane is consistent with the 
treatment in ref. [l3|. 

Overall, the observed temperature dependence is well 
reproduced by our theory down close to Neel temperature 
T/v ~ 2.3 [K], where the interchain interactions should be 
important also for ESR. Discrepancies are however found 
in the high temperature regime (> 200K) for all the three 
directions of H in ref. 14], and in the low temperature 
regime (< 50K) for H \\ c in ref. 0|. The former could 
be due to the low signal-to-noise ratio at high tempera- 
tures. The latter data may be also questionable for the 
absence of the anomaly (presumably related to the higher 
dimensional effects) below lOK, which are observed in all 
the other data set. 

The estimated values of the exchange anisotropy, espe- 
cially J'c/J ^ —0.15 for the data from ref. |13|, seems too 
large compared to the typical value of a few percent for 
a Cu ion. This could be due to a possible overestimation 
of the g-shift in the analysis in ref. 1^ to extract the 
direction-dependent effective (7-factor from the direction- 
averaged data on the powder sample. 

In any case, the temperature dependence of the g-shift 
of LiCuV04 is better accounted by the present theory 
than by the classical one |^, especially at low tempera- 
tures. We hope that the present study stimulates further 
experimental studies on the S" = 1/2 antiferromagnetic 
chain for more precise comparisons with the theory. 
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FIG. 3: Normalized effective g-shift 5g/{goof) ~ (ffcff — 
Qao) / {goaf) where / is the direction-dependence factor 0, 
is plotted against the temperature T. Points are experimen- 
tal data for LiCuV04 Hal and while the curve is the 
theoretical result lllll and the classical results 1141 . Inset: 
the behavior at low-temperatures. 
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